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ABSTRACT 



Anomalous X-ray Pulsars (AXPs) are now established to exhibit significant 
X-ray variability and be prolific glitchers, with some glitches being accompanied 
by large radiative changes. An open issue is whether AXP glitches are gener- 
ically accompanied by radiative changes, relevant for understanding magnetar 
physical properties. Here we report on an analysis of archival X-ray data from 
the AXP IE 1841-045, obtained between 1993 and 2007. This AXP, located in 
the center of SNR Kes 73, has exhibited three glitches between 2002 and 2007, as 
determined by RXTE monitoring since 1999. We have searched for evidence of 
phase-averaged flux variability that could be present if glitches in AXPs are usu- 
ally accompanied by radiative changes. We find no evidence for glitch-correlated 
flux changes from this source, arguing that such behavior is not generic to AXPs. 

Subject headings: pulsars: individual (IE 1841—045) — X-rays: stars — stars: 
neutron 



1. Introduction 



X-ray variability in AXPs is now established as a characeteristic property of this source 
class. Their flu xes and spectra va ry on a variety of different time sc ales, fr om a f e w mi l- 
liseconds (e . g. iGavriil et al.l 120021 ) . to several years ( jDib et al.l 120071 ) . See iKaspil ( 120071 ): 



Mereghettil ( 120081 ) for recent reviews of AXP variability. 



AXPs are also now known to be prolific glitc hers, and i nclud e some of the most active 
glitchers known in the neutron-star population (jDib et al.l 120081 ). By 'glitch' we mean a 
sudden rotational spin-up, or similar timing anomaly. One of the largest glitches seen so far, 
from AXP IE 2259+586, was accompanied by major radiative changes, including bursts and 
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(Kaspi et al. 


2003; 


Woods et al. 


2004) 



Dib et al.l (120091 ) showed that glitches in AXP IE 1048.1—5937 were also accompanied with 
radiative enhancements. Such events are thought to be the result of sudden yielding of the 
neutron-star crust due to internal stresses caused by the decay of the magnetar-strength 
field. The restructuring results in changes to the stellar interior - as evidenced by the glitch 
- and to the stellar exterior - as evidenced by the dramatic radiative changes. Observations 
of such AXP outbursts are a potentially powerful probe of the physics of magnetars (e.g. 



Eichler fc Shaisultanovll2010l ). 



Links between t he X-ray variabil i ty and glitches of AXP RXS J170849.0-400910 have 
also been reported ( IRea et al.l l2005t ICampana et al.l l2007t llsrael et al.l 120071 : iGotz et al. 
20071 ). These authors suggest the exi stence of a gener al correlation between magnetars' 
flux and glitch epochs. Specifically, IGotz et al.l (120071 ) reported ~40% flux changes in 
RXS J170849. 0—400910 based on eight observations made over the course of ~9 yr, dur- 
ing which they report four glitches. If correct, this suggests that glitches are usually, and 
possibly always, accompanied with radiative changes. However, the sparsity of observing 
epochs compared to the number of glitches for this AXP thus far is problematic in proving 
the variability is glitch-correlated. Additionally, no comparable pulsed flux changes were ob- 
served in the sa me time span, du ring which such measurements were available regularly on 
a monthly basis ( jDib et al.ll2008l ). This apparent conflict could, however, be explained if the 
pulsed fraction were precisely anti-correlated with flux. An anti-correlation between pulsed 



fracti on and flux has been seen in AXP IE 1048.1—5937 ( ITiengo et al.l l2005t iTam et al. 



20081 ). although not to a degree that render pulsed flux variations absent. 



Here we investigate the hypothesis that AXP glitches are generically accompanied by 
radiative changes by considering AXP IE 1841— 045 in supernova remnant (SNR) Kes 73, 
one of the most frequent glitchers among AXPs. From RX TE monitoring, we know that 
IE 1841-045 has had three glitches between 1999 and 2008 (bib et all 120081 ) . The glitches 
occurred on 2002 July 9, 2003 December 24, and 2006 March 29, and had Av/v of 5.63 x 
10" 6 , 2.45 x 10~ 6 and 1.39 x 10" 7 , respectively, and were not accompanied by any X-ray 
pulsed flux changes. If, as seen in IE 2259+586 and IE 1048.1—5937 and reported for 
RXS J170849. 0—400910 , glitches are generically accompanied by radiative changes, and, if as 
for RXS J170849. 0—400910 , such radiative changes are not necessarily apparent in the pulsed 
flux data, it is possible that the phase-averaged flux (unavailable from RXTE monitoring) 
varies in concert wit h glitches in IE 1841 —045 . We investigate this possibility here. Also 
worth noting is that iGotthelf et al.l (119991 ) studied the timing of AXP IE 1841— 045 using 
archival GINGA, ASCA, RO "SAT and RXTE data taken between 1993 and 1999 and found 
no glitches with Avjv > 5 x 10~ 6 . 
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In this paper, we report on our analysis of archival X-ray data for IE 1841—045 collected 
by ASCA, Chandra, XMM and Suzaku during the past 17 years, including two observations 
that were made fortuitously very closely following glitches. We have looked for correlations 
between the AXP's flux variability and glitch epochs. In section [2] we describe the observa- 
tions and data reduction process. In section [3] we describe our spectral analysis and AXP 
flux extraction method. Our results and conclusions are discussed in section HJ 



2. Observations 



For this study, we searched online X-ray archives for all existing observations of the 
Kes 73 field. We found a total of eleven observations from four different focussing X-ray 
observatories (listed in Table 1). We do not include in our analysis the many RXTE ob- 
servations, as due to its no n-focussing natur e, these provide only pulsed flux measurements, 
and are already published (IDib et al.l 120081 ). Next we report on our analysis of the eleven 
focussing-telescope observations. 



2.1. ASCA Observations 



Seven ASCA (ITanaka et al.lll994f ) observations of the AXP IE 1841-045 and the SNR 
Kes 73 were taken between 1993 and 1999 (see Table CD) . Our analysis beg an with the screened 
data from the two Gas Imaging Spectrometers (GISs. [Burke et al.lll994l ). which were filtered 
with the standard revision 2 screening criteria. Given the angular resolution of the GISs, the 
SNR (~2') was unresolvable in the images, therefore the spectra we extracted from the GISs 
contain photons from both the AXP and the SNR. Using the f tool xselect0, we extracted 
spectra from source regions of 9'. 8 radius (a region large enough to encircle the extended 
emission from AXP IE 1841— 045 and Kes 73) and background spectra from regions of radius 
~5', away from the source region for all the ASCA GIS observations. The GIS spectra were 
then combined with Redistribution Matrices File (RMF) of the GISs and Auxiliary Response 
File (ARF) files generated using the ftool ascaarf , and grouped with a minimum of 25 
counts per bin. Finally the exposure of the grouped spectra were corrected for the deadtime 
effect using the ftool deadtime. In this study, we di d not include the spectra extracted 
from the two Solid-state Imaging Spectrometers (SISs, IOhashilll996l ). primarily because of 
the significantly fewer counts collected by these instruments. 



1 http://heasarc.gsfc.nasa.gov/ftools/xselect/ 
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2.2. Chandra Observation 



The AXP and SNR were observe d by the Chandra X- ray Observatory with the Advanced 
CCD Imaging Spectrometer (ACIS, iGarmire et al.ll2003l ) in timed exposure (TE) mode on 
2000 July 23 and in continuou s clocking (CC) m ode on 2000 July 29 (Tabled]). The data 
were analyzed and reported by Morii et al. ( 2003 ). Here we used Ciao version 3.4§. Because 
the spectrum of the pulsar in TE mode was heavily affected by pile-up, we did not use these 
data; instead we extracted the spectrum of AXP IE 1841— 045 from the CC mode data. In 
these data, the image of the pulsar and SNR were collapsed into one dimension. From the 
level 2 event list provided by the ChandraX-ray Center, we extracted the spectrum of the 
pulsar using a box-shaped region capturing a 2". 5 long segment along the one-dimensional 
image and centered on the pulsar. The background spectrum was extracted from two 7". 5 
long segments adjacent to the source region. The resulting source and background spectra 
were then combined with RMF and ARF files generated using the psextract command and 
grouped with a minimum of 25 counts per bin. 



2.3. XMM Observations 



The AXP was o bserved by the Newton X-ray Multi-Mirror Mission {XMM ) observatory 
Hansen et aD boOlh on 2002 Octobe r 5 and 2002 Octob er 7 (Table Q} with the European 
Photon Imaging Camera (EPIC) p n (Striider et al.ll200lh camera operating in large window 
mode and the EPIC MOS cameras ( Turner et al. in full window mode. For our analysis, 

we used the XMM Science Analysis System (SAS) version 8. 0.d and calibrations (updated 
2008 Oct 3). Given the pn and MOS cameras' angular resolution, IE 1841—045 can be 
resolved from Kes 73. For the two XMM observations, we used only the data from the EPIC 
pn camera to take advantage of its larger photon collecting area and to avoid cross-calibration 
issues between the pn camera and the mos cameras. We extracted the pulsar's spectrum 
from a circular region of radius 32". 5 (a radius large enough to capture more than 90% of 
the photon events from the point source) centered on the pulsar. Background spectra were 
extracted from an annular region of radius between 35" and 115" centered on the pulsar. 
The pulsar spectrum was grouped with a minimum of 25 counts per bin and then combined 
with the background spectrum and RMF and ARF files generated by the SAS software. 



2 http://131.142.185.90/ciao3.4/index.html 
3 See http://xmm.esac.esa.int/sas/8.0.0/ 
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2.4. Suzaku Observations 



AXP IE 1841— 045 and Kes 73 were also observed by the Suzaku o bservatory (IMitsuda 



20071 ) on 2006 April 19 (Tabled]). The data analysis was reported by iMorii et al 



f l2008f ). 



Onboard Suza ku , there are two X-ray detectors: the X-ray Imaging Spectrometer (XIS, 
Koyamal 120071 . c onsisting of four CCDs sensitive in soft-X-ray band) and the Hard X-ray 
Detector (HXD. iTakahashil 120071 . sensitive to 10-600 keV X-rays). Here we present a spec- 
tral analysis of the XIS data only. Given the angular resolution of the XIS, the SNR was 
unresolvable in the XIS image. Therefore the spectra we extracted from the XIS detectors 
contain photons from both the AXP and the SNR. We used cleaned events screened by the 
standard pipeline processing version 2.0.6.130. The source spectra were extracted from a 
circular region of 260" radius. Background spectra were extracted from an annulus region of 
radius between 260" and 520". The extracted spectra were grouped with a minimum of 25 
counts per bin, and then combined with the RMF and ARF files generated using the f tools 
xisrmfgen and xissimarf gen. 



3. Spectroscopy 

Among the eleven observations we used, the AXP can be resolved from the SNR Kes 
73 only in the Chandra and XMM observations. For these, it is possible to extract either 
the neutron star's spectrum or the combined spectrum of the neutron star and the SNR. By 
contrast, only combined spectra can be extracted from the ASCA and Suzaku observations. 

The spectra of AXPs are often parametrized by a blackbody plus a power-law (although 
this is known to be an approximation to a likely Comptonized blackbody spectrum - see 
Thompson et al. 2002 and, for example, Rea et al. 2008). The spectra of SNRs are 
often fit with models like, for example, VSEDOV (a plane-parallel shock radiation model 
with separate ion and electron temperatures), VNEI (a non-equilibrium ioni zation collisional 



plasm a model), or VPSHOCK (a plane-parallel shocked plasma model). See lBorkowski et al. 



20011 for a review of these models. In this paper, we modeled the neutron star radiation 
with a blackbody plus power-law (BB+POW), and the SNR radiation with a VSEDOV 
model, using xspec 12.5.(0 The focus of our investigation is on the AXP's X-ray flux; 
in modelling the SNR we sought only a suitable parameterization to allow us to subtract 
off its flux reliably. As we show below, the VSEDOV model is adequate for these purposes. 



4 http: / /heasarc. nasa.gov / docs / suzaku/processing/ criteria_xis.html 
5 http: / /heasarc. nasa.gov /docs/ xanadu /xspec / 
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We modeled the interstellar absorption by multiplying a WABS (a photo-electric absorption 
model, in which the interstellar absorption is characterized by a single parameter Nh, the 
neutral hydrogen column density along the line of sight) model to both the BB+POW and 
VSEDOV models. 

Figure [I] shows the combined spectra from XMM and the components of the best-fit 
model. The AXP power-law component clearly dominates the spectra above ~4 keV. There- 
fore, we chose to study the AXP flux in the 4-10 keV band only, in order to minimize SNR 
contamination. Nevertheless, we still attempted to remove the remaining small contribution 
of SNR in this band for the fluxes measured from ASCA and Suzaku observations, so that we 
could compare them with the neutron-star-only fluxes measured with Chandra and XMM. 

It is reasonable to assume that both the interstellar absorption and SNR radiation do 
not change over a time scale of about a decade. Consequently, in our attempt to remove 
the SNR flux, we used the same Nh and VSEDOV parameters for all the spectra of the 
different observations, and allowed only the BB+POW model to vary from observation to 
observation. 

When fitting a WABS (BB+POW+ VSEDOV) model to those spectra containing emis- 
sion from both the AXP and SNR, it is challenging to constrain the normalization pa- 
rameters of both the BB and VSEDOV models simultaneously, because these two models 
dominate the same energy band and their parameters are highly covariant. Fortunately, 
the neutron star can be spatially resolved out in the XMM observations, so we can use 
them to determine the relative strength of the two spectral components. We therefore fit a 
WABS (BB+POW+ VSEDOV) model to the combined (AXP+SNR) XMM spectrum. To 
ensure that we had the correct BB+POW model for the neutron star, we simultaneously fit 
the spectrum extracted from only the neutron star, requiring common neutron-star spectral 
parameters. Hence we could determine the parameters and normalization of the VSEDOV 
model for the SNR. 

To further improve the spectral model, next we included the ASCA and Suzaku spectra 
and the Chandra CC-model neutron-star spectrum and performed a large joint fit, which 
required multiple iterations in order to converge. The result was a value for Nh and for 
the VSEDOV model parameters that fit all the SNR spectra reasonably well, albeit not 
perfectly (reduced \ 2 — 1-36 for 3801 degrees of freedom for the joint fit). The best-fit Nh 
and VSEDOV parameters can be found in Table [2j 

Note that for all the spectral fitting described above, we used the 0.8-10 keV band. 
However for Suzaku , we found that there were always significant residuals in the range 1.7- 
3.5 keV and above 9 keV, and in general, these residuals differed significantly among the 
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four XIS instruments. Therefore, we ignored the 1.7-3.5 keV band and above 9 keV for 
the Suzaku spectra. Furthermore, the SNR Kes 73 is larger than the field-of-view of the 
Suzaku XISs and was not entirely captured. Therefore, when fitting the Suzaku spectrum, 
we allowed the VSEDOV normalization parameter to vary. 

Finally, by using the best-fit Nh and VSEDOV model, we could remove the SNR flux 
contribution from the ASCA and Suzaku observations in separate fits to their spectra, hence 
measuring the AXP fluxes. For the Chandra and XMM observations, we simply fitted the 
resolved AXP spectra with an absorbed BB+POW model. The best-fit BB+POW param- 
eters and the measured 4-10 keV unabsorbed fluxes of the AXP are presented in Table [3j 
The fluxes are also plotted in Figure [2j The uncertainties we report on the 4-10 keV un- 
absorbed AXP fluxes were estimated by the measured fractional uncertainties on the 4-10 
keV absorbed total fluxes. For those observations in which the best-fit reduced x 2 was larger 
than unity, we multiplied the reported flux uncertainties by the square root of the reduced 
X 2 ■ This is to account for systematic errors in our imperfect modelling of the SNR spectrum. 

To confirm the robustness of the measured AXP flux values, we repeated the entire above 
analysis by modelling the SNR spectra in different ways. We found that both the VNEI and 
VPSHOCK models can fit the SNR spectra as well as the VSEDOV model, and that the 4-10 
keV unabsorbed neutron star fluxes we measured using these models are consistent within 
reported uncertainties with the values we found using the VSEDOV model. Therefore, we 
feel confident that the phase-averaged AXP flux values we report are well constrained and 
robust. However, the BB+POW spectral parameters we measured when fitting the SNR 
with different models were not as robust as the 4-10 keV fluxes, changing significantly with 
SNR model. As we do not believe them to be reliably determined, we do not quote their 
uncertainties in Table [3j and do not consider them further. 

Also in Figure [2J we show the AXP's 2-10 keV pulsed flux as measured in monitor- 
ing observations with RXT E since early 199 9. Details about how these pulsed fluxes were 



determined are provided in iDib et al.l (120081 ). The RXTE pulsed fluxes show no significant 



variations. The vertical lines in Figure [2] indicate the epochs of the three observed glitches 



(IDib et alJl2008h . 



4. Discussion and Summary 

The goal of this study was to see whether the prolific glitching AXP IE 1841—045 shows 
phase-averaged flux variability, in spite of showing no evidence for pulsed flux variability. 
Also, we wished to determine whether any variability is correlated with its glitches as has 
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been seen in AXP IE 2259+586 in its 2002 major outburst, in IE 1048.1-5937, and also 
reported for RXS J170849.0-400910 . 

As is clear from Panel a of Figure |21 in the 4-10 keV band, the neutron star's flux 
did not vary by more than ~30% in 13 years. Interestingly the largest variations we find 
are in the multiple pre-1999 ASCA observations: in those seven observations, a fit to a 
constant flux results in a reduced x 2 of 3.7 for 6 degrees of freedom, which has a probability 
of occurring by chance of — 0.001. However, duri ng this time, there were certainly no large 
glitches (Az//z/ < 5 x 10~ 6 ) (fcotthelf et al.l Il999h . 



The ASCA fluxes all appear to be higher than the fluxes measured from other observa- 
tions. If this difference is accurate, then the AXP's flux must have dropped around 2000, 
and we might in principle expect some change in the pulsed flux as monitored by RXTE at 
that epoch. However, the pulsed flux of the pulsar was constant around that time as seen in 
Panel b of Figure [2j Therefore, we suspect that the relative increase in the flux in the ASCA 
observations compared with those of the other observatories may be due to instrumental 
calibration issues. 

The fluxes measured from the last four observations taken by Chandra , XMM and 
Suzaku can be fitted with a constant flux model (reduced x 2 — 1-9 f° r 3 degrees of freedom, 
corresponding to a probability of having occurred by chance of 0.125). Thus we conclude 
that the phase-averaged 4-10 keV fluxes of the last four observations were consistent with 
being constant, and we put an upper limit of 11% on long-term variability in this energy 
band. Importantly, the two XMM observations were taken only 88 and 90 days after the first 
glitch, and the Suzaku observation was taken only 27 days after the third glitch. By contrast, 
the 4- 10 keV flux of IE 2259+586 was 50% higher than in quiescence 21 days after its 2002 
glitch rtzhu et al.ll2008h. t hat of IE 1048. 1-5937 was a factor of 6 higher >38 days after its 
2007 glitches JTam et al.l l2008h . and was 50-7 0% higher f o r RX S J170849.0-400910 -53 



days after its first 2005 glitch as inferred from iGotz et al. Therefore, we conclude 

that unlike in IE 2259+586, IE 1048.1-5937 and possibly RXS J170849.0-400910 , there is 
no evidence for glitch-correlated flux changes in AXP IE 1841—045 . 

One caveat of our study is that we were limited to the harder part of the neutron 
star's emission spectrum. The flux from the blackbody component was not well constrained. 
Therefore, we cannot rule out changes in the neutron star's thermal radiation, only changes 



6 This is calculated based on the reported 1-10 keV fluxes of RXS J170849.0-400910 from its 
2003 XMM observation and 2005 Swift observation. We assumed a power-law to blackbody flux 
ratio of 3 for the 1-10 keV fluxes, and then calculated the 4-10 keV fluxes using webPIMMS 
(http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html) 
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in the power-law component in the 4-10 keV band, which we note constitutes 
stellar flux (BB+POW) in the 1-10 keV band. 



-0.25 of the 



Thus our results argue against the hypothesis that there exists a generic correlation 
between the X-ray flux variability and glitch epochs in AXPs, further supporting the argu- 



ment that glitches in AXPs can be either radiatively loud or radiatively silent (IDib et al. 



20081 ). There is of course precedent for radiatively silent glitches in neutron stars, in that 
no rotation-powered pulsar glitch has ever been reported to be accompanied with any ra- 
diative change, alth ough rapid X-ray follow-up has only been accomplished in one case 
( IHelfand et al.ll200ll ). Any physical model of magnetar glitches will have to explain the si- 
multaneous existence of both types. This is true of even a single source, as there is evidence 
that A XP IE 2259+586h as both, gi yen that its most recen t glitc h showed no pulsed flux 
change ( IDib et al.ll2008bl ). Recently lEichler fc Shaisultanovl (120 10f ) have argued that AXP 
glitches are triggered by energy releases at depths below ~100 m in the crust, with angular 
momentum vortex unpinning being due to global mechanical motion triggered by the en- 
ergy rele ase, not by heat as has been proposed in the context of rotation-powered pulsars 
glitches ( iLink fc Epstein! Il996l ; iLink fc Cutlerl 120021 ) . If mechanical triggering occurs, then 
radiatively silent glitches of the same amplitude as radiatively loud glitches are possible since 
less ener gy is required to trigger a glit ch than to cause a substantial X-ray brightening. If 
so, then lEichler fc Shaisultanovl (120101 ) predict that all AXP glitches should occur simulta- 



neously with or before the observed X-ray brightening; this can be tested with continuous 
(daily or better) X-ray monitoring observations. Moreover, although mechanical unpinning 
of vortices by activity in the lower crust does result in heat release, the latter could take 
as much as several years to reach the surface. This could help explain the lon g-term X-ray 
variability trends that have been reported in some AXPs (e.g. IDib et al.l 120071 ). 
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Table 1. X-ray observations of AXP IE 1841— 045 used in this study. 



Date 


Observatory 


t a (ks) 


Offsct b (') 


1993 Oct 11 


A SCA 


40 


6.8 


1997 Apr 21 


ASCA 


9 


7.3 


1998 Mar 27 


ASCA 


39 


6.5 


1999 Mar 22 


ASCA 


20 


5.3 


1999 Mar 29 


ASCA 


20 


5.2 


1999 Apr 06 


ASCA 


19 


5.2 


1999 Apr 13 


ASCA 


21 


5.2 


2000 Jul 29 


Chandra 


10 


0.097 


2002 Oct 05 


XMM 


3.8 


1.152 


2002 Oct 07 


XMM 


4.4 


1.144 


2006 Apr 19 


Suzaku 


98 


3.9 



a The effective exposure time of the instrument used 
for the spectral analysis in this paper. 

b The pointing offsets of the observations relative to 
the position of the AXP. 
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Table 2. WABS and VSEDOV parameter values 



Parameter 


Value 


N H 


2.77 xl0 22 cm- 2 


kT a 


0.56 keV 


kT b 


0.56 keV 


Mg 


0.63 a 


Si 


1.10 


S 


1.92 


Ca 


0.32 


Fe 


0.57 


Ni 


2.51 


T 


1.50 x 10 11 s 


Redshift 


0.00 


Norm b 


0.33 


H...Ar c 


1.0 



a The clement abundances 
quoted here are the relative 
abundances based on the Solar 
mixture abundances. 

b Normalization parame- 
ter for the VSEDOV model, 

■■ rn rr, M2 r where Da 

is the angular diameter distance 
to the source (cm), and n ei nn 
(cm -3 ) are the electron and 
hydrogen densities respectively. 

c Elements H...Ar: H He C N 
O Ne and Ar were fixed to the so- 
lar abundance because the SNR 
spectra are not sensitive to them. 
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Table 3. Measured AXP model parameters and fluxes. 



Date 


Observatory 


MJD 


kT a (keV) 


pa 


F NS h (10- n erg cm^s- 1 ) 


1993 Oct 12 


ASCA 


49272.1 


0.20 


2.72 


1.16(4) 


1997 Apr 21 


ASCA 


50559.6 


0.19 


3.00 


1.2(1) 


1998 Mar 27 


ASCA 


50899.95 


0.20 


2.70 


1.36(5) 


1999 Mar 22 


ASCA 


51259.7 


0.20 


2.89 


1.29(4) 


1999 Mar 30 


ASCA 


51267.0 


0.19 


2.80 


1.39(5) 


1999 Apr 06 


ASCA 


51274.6 


0.21 


2.74 


1.39(5) 


1999 Apr 14 


ASCA 


51282.3 


0.21 


2.68 


1.36(4) 


2000 Jul 29 


Chandra 


51754.32 


0.40 


2.10 


0.975(3) 


2002 Oct 05 


XMM 


52552.16 


0.22 


2.59 


1.05(4) 


2002 Oct 07 


XMM 


52554.16 


0.21 


2.47 


1.08(3) 


2006 Apr 19 


Suzaku 


53844.8 


0.38 


1.94 


1.03(2) 



a The reported best-fit kT and T parameters of the BB+POW model vary depending on 
the assumed SNR model; we do not report the fit uncertainties of these parameters as they 
do not reflect the true uncertainties. The values, determined while assuming the VSEDOV 
model for SNR, are provided for reference only. 

b Phase-averaged unabsorbed neutron-star flux in the 4-10 keV band. The number in 
parenthesis represents the la uncertainties in the last digit. See text for details. 
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Fig. 1.— The unfolded XMM pn spectra and components of the best-fit 
BB+POW+VSEDOV model. 
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Fig. 2.— (a) 4-10 keV unabsorbed phase-averaged fluxes from AXP IE 1841—045, as 
determined by our analysis, (b) 2-10 keV pulsed fluxes from the pulsar as me asured by 



RXT E , in units of counts per second per Proportional Counter Unit (PCU) (see iDib et al. 



20081 . for details). 



